Physical Chemistry Chemical Physics

The process of drug unbinding is of immense importance in the field of biophysics and therapeutics. The behavior of these systems is greatly influenced by their thermodynamic and kinetic properties. Therefore, it is crucial to accurately estimate the ligand binding free energies and rate of ligand dissociation, yet these processes are often governed by rare event transitions that lie beyond the reach of standard brute-force molecular dynamics simulations. While enhanced sampling simulations offer a solution, their efficacy is strictly contingent upon the selection of appropriate collective variables (CVs) which is non-trivial for complex systems like protein-ligand complexes. In this study, we present a method to derive optimized CV from transition state region (TS) via an interpretable machine learning (ML) model, Elastic Net. By employing some physically intuitive order parameters, the derived optimized CV from the TS-region greatly accelerate ligand binding-unbinding transitions and achieves rapid free energy surface (FES) convergence across diverse systems including buried and solvent exposed active sites such as Trpsin-benzamidine complex, host-guest systems and sodium epoxidase etc. Intriguingly significant contribution of the ligand hydration is found in the optimized CV which depicts crucial role of solvent in driving ligand binding-unbinding transitions. The estimated binding free energies for different protein-ligand complexes match quite well with experiments, while maintaining a low computational cost. The derived optimized CV is also used to calculate the ligand residence times across different systems and calculated residence times are within the experimental range for all systems, again with very little computational costs. Moreover, we show that the optimized CV constructed from TS region via an interpretable ML model is transferable across diverse systems, offering a robust and scalable framework for drug discovery and investigation of complex biomolecular recognition.

The High Mobility Group Box 1 (HMGB1) protein performs multiple essential functions in the body, ranging from DNA regulation to the activation and mediation of immune responses. However, HMGB1 has been also implicated in several pathological conditions, such as rheumatoid arthritis, sepsis, autoimmune diseases, tumors, and Alzheimer's disease. As a result, HMGB1 is of increasing interest as a therapeutic target. Binding to heparin has been reported to inhibit HMGB1's pathological activity during sepsis in clinical settings. In this work, we compare the interactions of HMGB1 with heparin and its' synthetic analog linear polyglycerol sulfate (lPGS) from the viewpoint of stability and changes to association behavior. This analysis focuses on thermal stability, secondary-structure changes, and particle-size evolution using nano-differential scanning fluorimetry (nanoDSF), circular dichroism spectroscopy (CD), and dynamic light scattering (DLS).

1.Surface functionality and pH play a decisive role in governing the structural dynamics, hydration, and drug-binding behaviour of dendrimers. Here, all-atom molecular dynamics (MD) simulations were performed on five generations of PAMAM (G1-G5) and PETIM (G2-G6) dendrimers with O-core and N-core architectures, functionalized with amine, carboxylic acid, or sugar terminal groups under different protonation states. Protonation of the tertiary branch-point amines expands the dendrimer structure, increases internal porosity and hydration, and enhances structural fluctuations across both families. In contrast, non-protonated amine -NH2 (NP) and carboxylic acid -COOH (NP) terminated dendrimers, together with deprotonated carboxylate-COO- (DeP) systems, retain comparatively compact conformations. Sugar-functionalized dendrimers ({beta}-galactose-terminated PETIM and D-glucose-terminated PAMAM) are most hydrated and structurally rigid, whereas amine-terminated dendrimers exhibit the greatest conformational dynamics. PAMAM dendrimers with -NH2, -NH3+, and -COO- terminal groups are generally more hydrated than their PETIM counterparts. However, {beta}-galactose-terminated PETIM dendrimers are more hydrophilic than D-glucose-terminated PAMAM dendrimers. N-core PETIM dendrimers also adopt more compact and spherical conformations than equivalent O-core PETIM dendrimers. Drug-binding MD simulations show that curcumin binding is dominated by van der Waals (vdW) interactions, whereas doxorubicin complexation is primarily driven by electrostatic interactions. Among the investigated surface functionalities, -NH2 (NP), -NH3+ (P), -COOH (NP), and -COO- (DeP) terminations exhibit the most favourable drug-binding characteristics. Except for deprotonated carboxylate systems, curcumin binds more strongly than doxorubicin. Overall, these findings establish molecular-level relationships between surface functionality, protonation state, dendrimer architecture, and drug-binding behaviour, providing design principles for pH-responsive dendrimer nanocarriers with enhanced drug-loading and controlled-release performance. TOC O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=75 SRC="FIGDIR/small/742721v1_ufig1.gif" ALT="Figure 1"> View larger version (35K): org.highwire.dtl.DTLVardef@119bf29org.highwire.dtl.DTLVardef@1554d86org.highwire.dtl.DTLVardef@154a254org.highwire.dtl.DTLVardef@16d5c5b_HPS_FORMAT_FIGEXP M_FIG C_FIG

Photoprotection against excess energy is essential for the survival of photosynthetic organisms under adverse conditions. In plants, excess energy can be dissipated as heat through non-photochemical quenching (NPQ) of chlorophyll fluorescence, involving the trimeric light-harvesting complex II (LHCII), the major antenna of photosystem II. How NPQ affects antenna proteins remains debated, especially as most studies focus on short-lived components artificially induced in vitro. Here, we characterize the effects of qH, a long-lived NPQ component, on the fluorescence properties of natively quenched LHCII. Single-molecule fluorescence measurements, combined with biochemical and biophysical ensemble approaches, reveal a larger and more quenched subpopulation of LHCII trimers exhibiting fluorescence intermittency in samples with qH compared to those without. This behavior is linked to a small conformational change that stabilizes a quenched state, enhancing photoprotection at the antenna level. These findings provide new insights into sustained NPQ and its role in regulating energy dissipation under natural light conditions.

This present work shows that E-hook fragments possess functional structure differences governed by electrostatic interactions and sequence composition. The acidic C-terminal tails of tubulin, known as E-hooks, play a central role in regulating interactions between microtubules and motor proteins, microtubule-associated proteins, and enzymatic modifiers. Despite their functional importance, the intrinsic structural properties of these peptide segments remain poorly characterized due to their intrinsically disordered nature. In this work, we present quantum-mechanically optimized structures of hexamer peptides derived from {beta}-tubulin E-hook sequences. Density functional theory calculations were used to optimize peptide geometries using progressively larger basis sets. From the optimized geometries we calculated theoretical Raman spectra, Ramachandran backbone dihedral distributions, and measured radii of gyration to resolve composition dependent structural tendencies. The combined Raman and conformational analyses provide a systematic computational approach for comparing simulated and experimental Raman spectra of tubulin E-hooks and other intrinsically disordered proteins and offer insight into how E-hooks contribute to the recognition mechanisms underlying the tubulin code.

The rising incidence of invasive fungal diseases emphasizes the need for novel therapeutic strategies, including membrane-targeting antifungal agents, which require representative lipid models for detailed molecular-level studies. In this work, we propose a consensus quinary fungal plasma membrane model based on lipidomic literature data, specifically PC:PE:PI:PA:PS phospholipid model with ratio of 44:29:13:8:6. Using a bottom-up approach, we characterized the biophysical properties of this system - with particular emphasis on mechanical parameters such as bending rigidity and area compressibility - by combining molecular dynamics simulations with experimental flicker-noise and ATR-FTIR spectroscopies. Furthermore, we investigated the effect of two key non-phospholipid components: ergosterol and triacylglycerols. Biophysical analysis revealed that DPPI and its specific interactions with DSPS induced the most substantial deviations in baseline membrane parameters, particularly area per lipid, membrane thickness, and area compressibility, while DSPS influenced bending rigidity change and DLiPA primarily affected lipid packing defects. In addition, ergosterol and TGs were found to influence all of the investigated parameters to different degree. Notably, the overall biophysical profile of the proposed FPMM closely mimicked that of natural vesicles derived from yeast lipid extracts, establishing this model may provide a reliable platform for studying fungal membrane biophysics and lipid-targeting interactions.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

Nanodiscs and bicelles are widely used as membrane mimetics for structural studies of membrane-associated systems. Studies have reported that their magnetic alignment behavior and phase stability are highly sensitive to composition and temperature. In this study, we systematically investigate the effects of cholesterol on bicelle formation and magnetic alignment in DMPC + 0.2 glycyrrhizic acid (GA) systems using a combined 31P and 14N solid-state NMR experimental and simulation-based approach. Temperature dependent 31P NMR spectra reveal a clear transition from vesicle dominant to aligned bicelles/nanodsics phase, while 1N quadrupolar splitting and lineshape analysis provides quantitative insights into heterogeneous lipid bilayer populations, distinguishing large aligned nanodiscs (B(L)), small nanodiscs (B(S)), and isotropic/random components (B(R)). A strong correlation is observed between the 31P derived bicelle fraction and the 14N B(L) population, confirming that macroscopic alignment in the presence of an external magnetic field directly reflects the growth of large, well-ordered nanodiscs. Cholesterol is found to play a critical dual role by modulating membrane order and curvature. At low cholesterol concentration (0 to 5 mole percent), nanodiscs alignment occurs gradually with increasing temperature, while at higher cholesterol concentration (15 to 25 mole percent), the alignment is delayed and accompanied by broader spectral features, indicating structural heterogeneity. Notably, 10 mole percent cholesterol consistently provides the optimal balance, enabling efficient temperature dependent conversion to aligned bicelles while maintaining high B(L) populations (about 70-80 percentage) and minimal isotropic fractions. In contrast, higher cholesterol maintains significant B(S) and B(R) populations, even at elevated temperature. The 14N quadrupolar coupling (Cq is approximately 8.5 to 9.2 kHz for aligned nanodiscs) remains nearly invariant across compositions, showing that cholesterol does not change local headgroup dynamics but instead redistributes lipid populations. These findings establish a combined 31P and 14N solid -state NMR approach provides a valuable platform for quantitatively correlating membrane structure, dynamics, and alignment, offering practical guidelines for optimizing bicelle systems for high resolution solid-state NMR studies of membrane associated biomolecules.

N-linked glycosylation is known to modulate the catalytic function of human carbonic anhydrase (HCA) IX, yet its influence on the underlying free-energy landscape remains largely unexplored. In the present work, we combine extensive all-atom molecular dynamics simulations with kinetic transition network analysis to investigate the effect of glycosylation on the conformational organization of the catalytic domain of HCA IX in both monomeric and dimeric forms. The multidimensional conformational space is discretized into distinct free energy minima using the distribution of reciprocal interatomic distances (DRID), and the effective barriers separating them are estimated using the max flow-min cut formalism. The corresponding free energy landscapes are visualized in terms of disconnectivity graphs, which provide a faithful representation of underlying kinetics. Minimum free energy paths, mean first passage times, as well as frustration metrics are computed to further quantify the effect of glycosylation on landscape topography. Unglycosylated systems are found to exhibit predominantly funnel-like landscapes, with a limited number of metastable states in the vicinity of the native protein fold. In contrast, glycosylation enhances landscape complexity, resulting in a wide array of relaxation timescales. Strikingly, the two glycan chains affect the landscape topography in distinct ways, despite having closely matching sequences. Dimerization couples the glycan chain dynamics, with transitions between key metastable states involving coordinated motions of both the chains. Our work illustrates that interpretation in terms of disconnectivity graphs and transition networks could reveal important insights into the organization of glycoprotein energy landscapes.

Cardiolipin (CL) is a phospholipid found in the inner mitochondrial membrane (IMM) where it increases the efficiency of ATP regeneration. We have investigated the hypothesis that this increase may result in part from CL concentrating H+ at the IMM surface through electrostatic interactions as the CL polar head is a dianion at physiological pH. To this aim, we compared the concentrations and movements of H+ at the surface of giant planar phosphatidylcholine (PC) membranes and 20% CL enriched PC membranes by recording their surface pH with the membrane-grafted pH probe fluorescein DHPE. CL enrichment of the membranes increased their surface H+ activity by a ~4 factor. Moreover, we observed non-gaussian spatial H+ concentration profiles with distance from a point H+ source with both PC and CL membranes suggesting that both lipids also induce interactions between probe molecules. A whole bath pH variation revealed that these interactions allow the traveling of reversible acidification fronts with constant speed over the membrane between high and low pH states. A reaction-diffusion model of these observations suggests that membranes support these fronts through a mechanism of autocatalytic (de)protonation of the membrane surface. In mitochondria, these fronts would result in transitions between high and low pH states, the low one having a larger H+ concentration in CL-enriched regions of the IMM. Such an increase at the inner leaflet of the IMM may increase efficiency of the respiratory chain whereas the increase at the outer leaflet may boost the ATP synthase rate.

A comprehensive understanding of the mechanisms underlying biological systems requires correlative analysis of multiple complementary molecular properties through multidimensional measurements. High-speed atomic force microscopy (HS-AFM) is a powerful tool for elucidating biomolecular structural dynamics at the single-molecule level with high spatiotemporal resolution. However, structural information alone is often insufficient for fully understanding the biological function mechanisms. Here, we report high-speed atomic force-Raman microscopy (HS-AFRM), which enables multimodal measurements combining video-rate structural imaging with chemical-bond analysis. Raman spectroscopy is a powerful, non-invasive technique that probes molecular vibrations to provide chemical information. We achieved several key technical developments that facilitated the seamless integration of HS-AFM and micro-Raman spectroscopy, allowing reliable correlative measurements of structural and chemical information. We validated the versatility of the developed system using representative samples, including two-dimensional materials and a protein. Furthermore, we demonstrate probing of changes in the surrounding environment, which are inaccessible by HS-AFM alone. Multimodal measurements incorporating fluorescence spectroscopy were also demonstrated as an additional practical extension. This multimodal approach substantially enhances the analytical capability of HS-AFM, providing a powerful platform for revealing correlated structural and chemical properties across diverse research fields.

Within the secondary cell walls of vascular plants, cellulose, hemicelluloses, and lignin associate via various covalent and non-covalent linkages to form an intricate supramolecular assembly. Although the chemical structures and cross-linking levels of the lignin- hemicellulose matrix exhibit substantial diversity in planta, precisely how these structural variations affect lignocellulose supramolecular assembly and macroscopic biomass properties remains largely elusive. Here, we conducted a comparative multi-scale ultrastructural analysis across engineered rice lines with targeted modifications in lignin aromatic composition and ferulate (FA)-mediated cell wall cross-linking levels. Combined solid-state nuclear magnetic resonance and wide-angle X-ray diffraction analyses revealed that depleting FA cross-linking disrupts cellulose crystalline structure and accelerates molecular mobility markedly more severely than altering the guaiacyl-to-syringyl (G/S) lignin ratio, generating a more loosened lignocellulose network. Small-angle X-ray scattering analysis further demonstrated that specific FA-depleted lines, but none of those with an altered G/S ratio, also exhibited disruptions in the nano- to mesoscale organization of cellulose microfibrils. Furthermore, FA- depleted lines generally displayed greater improvements in saccharification efficiency and more rapid thermal softening than G/S-lignin-altered lines, suggesting that disruptions in lignocellulose molecular assembly induced by FA depletion can broadly translate into macroscopic biomass properties. These findings establish a molecular basis for the pivotal role of FA cross-linking in dictating grass cell wall architecture, offering a promising structural target for advancing grass biomass utility and crop design.

Calcium oxalates are naturally occurring minerals, with the monohydrate form CaC2O4xH2O (COM) being the most stable. COM contains two crystallographically different water sites. We investigate the details of water internal dynamics in the high and low temperature phases of COM using 17O central transition solid-state NMR line shapes, as well as laboratory and rotating frame relaxation rates. The measurements were performed either under static or magic angle spinning conditions and in a wide temperature range from 343 to 180 K. The combination of all measurements allows for precise constraints on motional mechanisms, rate constants, and amplitudes of motions. The high temperature phase is dominated by large-angle fluctuations with an amplitude of about 100 degrees, identical in both sites. During the phase transition between 323 to 300 K, these large-angle jumps freeze out in one of the water sites, while remaining active in the other. In the low temperature phase from 280 to 180 K, small-angle fluctuations of 2-8 degrees in amplitude dominate the relaxation. Transverse relaxation rates also point to the existence of a very slow collective rocking motion down to about 220-200 K.

Antimicrobial resistance motivates the development of approaches capable of probing nanoparticle-bacterium interactions with nanoscale sensitivity. Here, synchrotron infrared nano-spectroscopy (SINS) is applied to investigate interactions between carbohydrate-coated silica nanoparticles and the Gram-negative model bacterium Escherichia coli at the single-cell level. Silica nanoparticles (SiO2) were coated with mannose, maltose, or trehalose to evaluate how surface carbohydrate chemistry influences their interactions with the bacterial envelope. Correlative electron microscopy revealed pronounced association of carbohydrate-SiO2 with the bacterial envelope, with features consistent with localization within the periplasmic region, whereas bare-SiO2 showed no detectable association. SINS measurements acquired directly on bacterial cells and at bacterium-nanoparticle interfaces revealed distinct, carbohydrate-dependent spectral signatures. Quantitative analysis of the amide I band used the I/I{beta} ; ratio, which describes the relative contributions of -helical and {beta}-sheet protein secondary-structure components, together with interface-dependent band-position analysis to characterize local spectral perturbations. Carbohydrate-SiO2 produced systematic changes in the I/I{beta} ; ratio, including at locations where nanoparticles were not directly observed, indicating that their effects extend beyond the sites of nanoparticle association. Comparison of measurements acquired on bacterial surfaces and at bacterium-nanoparticle interfaces further revealed that carbohydrate chemistry modulates both the magnitude and spatial extent of these spectral perturbations. Trehalose-SiO2 produced the largest interface-dependent amide I band shifts and a spectral component consistent with random-coil structures. Overall, these results demonstrate that carbohydrate surface chemistry modulates nanoscale protein conformational perturbations at the nano-bio interface and highlight SINS as a powerful approach for resolving chemically localized molecular responses at single-cell interfaces.

The recent development of Volumetric Imaging via Photochemical Sectioning (VIPS) has enabled nanoscale imaging of whole-mount tissue samples of virtually any size by embedding intact tissue in a photocleavable, superabsorbent hydrogel. However, the efficacy of sample embedding, imaging, and photochemical sectioning is fundamentally governed by the mechanical stiffness, structural stability, and photodegradation kinetics of the photocleavable hydrogel (PC-gel) polymer network. To elucidate the effect of the photosensitive crosslinker design on these critical properties, we synthesized a set of photocleavable crosslinkers (PCs) with varying polyethylene glycol (PEG) backbone lengths and prepared the corresponding PC-gels under a fixed monomer formulation and polymerization condition. We quantified and compared the viscoelastic properties of the formed PC-gels at their swollen states, and found that the crosslinker length markedly reshaped the PC-gel mechanics. In addition, we evaluated the light-triggered degradation of the PC-gels using both wide-field and spatially-controlled illumination. We found that PC-1000, PC-1500, and PC-2000 gels remained comparably photodegradable, all enabling on-demand, spatially confined decrosslinking under such illuminations. These results provide practical guidelines for modulating the crosslinker architecture of PC-gel polymer networks to achieve optimal physicochemical properties for whole-mount tissue imaging using VIPS.

Copper is an essential redox-active micronutrient, but agricultural soils are increasingly contaminated by copper from mining, industrial discharge, and intensive agrochemical use, pushing concentrations beyond levels plants can tolerate. Excess copper triggers Fenton-like reactive oxygen species (ROS) generation, mitochondrial dysfunction, and impaired growth. Existing mitigation strategies, such as soil amendments, phytoremediation, antioxidants, and different chelators, have been explored to reduce copper toxicity, but their effectiveness can be limited by immobilization, poor specificity, and environmental persistence. The present work introduces a nanoparticle-based strategy for the direct sequestration of excess copper coupled with protection against the oxidative damage caused by copper stress. Here, we report MPA-iron oxide nanoparticles (MIONPs), sequentially functionalized with chitosan, glutathione, and 3-mercaptopropionic acid, designed to simultaneously scavenge ROS, restore redox homeostasis, and chelate copper via surface thiol groups. MIONPs showed a significant increase in copper binding capacity over bare iron oxide nanoparticles (BIONPs) and, in copper-stressed Solanum lycopersicum seedlings, significantly improved germination and root/shoot growth, reduced intracellular ROS, restored mitochondrial membrane potential, and preserved nuclear integrity. This integrated design establishes MIONPs as a promising, dual-function nanoplatform for sustainable copper stress management in agriculture.

Photosystem II (PSII) is in eukaryotic phototrophs is generally considered to operate within a more restricted spectral range than photosystem I (PSI), in which long-wavelength chlorophylls are a well-established feature of the peripheral antenna. Whether eukaryotic PSII can acquire comparable far-red-associated properties through lineage-specific antenna diversification has remained unclear. Here we present a 3.09 [A] cryo-electron microscopy structure of the C2S2M2L2 PSII supercomplex from Euglena gracilis, a euglenophyte species harbouring a secondary plastid and unusual light-harvesting system. We show that the euglenophyte-specific antenna protein LhcE9 occupies the position corresponding to canonical Lhcb5, but in a markedly different orientation that creates a distinct interface with the PSII core, particularly with CP43. Combined structural, spectroscopic, mutagenesis and proteomic analyses support LhcE9 as the stably bound PSII antenna subunit most closely associated with the far-red state in the supercomplex. Excitation-energy-transfer calculations further indicate two fast lineage-specific antenna-to-core routes mediated by LhcE9 and PsbX. Together, these findings reveal an unexpected mode of PSII antenna diversification and provide a structural framework for far-red-associated light harvesting in PSII.

Residual dipolar couplings (RDCs) obtained via magnetic field-induced alignment offer a powerful, media-free approach for the structural analysis of biomolecules. However, their detection in short, fast-tumbling nucleic acids remains elusive at conventional magnetic fields due to insufficient alignment and sensitivity. Here, we demonstrate the direct observation of these RDCs at 1.3 GHz in a 14-mer hairpin fragment derived from an HIV-1 Vif-targeting aptamer. The 1JNH scalar couplings of imino protons were measured at fields ranging from 600 MHz to 1.3 GHz. While the coupling constants remained invariant between 600 and 900 MHz, a clear deviation was exclusively captured at 1.3 GHz for all base-paired stem residues, demonstrating the first media-free detection of field-induced RDCs in a short RNA of this size. This breakthrough arises from a synergistic B07/2 scaling, combining enhanced alignment ({propto} B02) and sensitivity ({propto} B03/2). These RDCs showed excellent agreement with the NOE-derived structure. Additionally, the flexible loop residue G8 exhibited no detectable RDC, but displayed a field-dependent TROSY/anti-TROSY intensity inversion at 1.3 GHz, reflecting an unusual 1H chemical shift anisotropy (CSA) tensor that corroborates the local base-packing environment. Our findings highlight 1.3 GHz NMR as an indispensable tool for the structural analysis of short RNAs.

Weak, reversible interactions underpin biomolecular recognition, and single-molecule fluorescence (smF) imaging techniques can provide unprecedented insight into those biological processes. Unfortunately, such studies often require micromolar concentrations of fluorophore-labeled biomolecules, which is beyond the accessible range of conventional smF microscopies. Here, we describe a surface-functionalization method based on cloud-point polyethylene glycol (PEG) grafting that enables widefield smF microscopy measurements at micromolar concentrations without the use of nanophotonic devices. Using conventional total internal reflection fluorescence (TIRF) microscopy, we detected single-molecule fluorescence resonance energy transfer (smFRET) from surface-tethered, donor-labeled target molecules with up to 8 micromolar concentrations of freely diffusing, acceptor-labeled analyte molecules in the background--two orders of magnitude higher than typical studies in the literature. Weak, DNA-hybridization and protein-RNA binding equilibria were measured across micromolar range titrations. Together with advances in high-background data analysis, the robust method presented here enables kinetic and thermodynamic analyses of weak biomolecular interactions, especially those limited by nonspecific adsorption and high fluorescence backgrounds, using only standard smF instrumentation.

The sequence-dependent mechanical landscapes of double-stranded nucleic acid (dsNA) remain largely unexplored beyond canonical dsDNA. We describe cgNA+, a coarse-grained predictive model of the mechanics of dsRNA, DNA:RNA hybrids, and epigenetically modified dsDNA, all parameterised from 1.26 milliseconds of atomistic simulations. cgNA+ predicts non-local sequence-dependent equilibrium shape and stiffness with errors an order of magnitude smaller than sequence-variability, while enabling exploration of numbers of sequences inaccessible to atomistic simulation. We show that dsNA equilibrium shape is strongly influenced by flanking sequence up to octamer context, with flexible dimer-steps more context-sensitive. CpG-modification alters equilibrium shape comparable to changes caused by single-nucleotide polymorphisms. Groove width analysis across dsNA decamers reveals strong sequence dependence, reflecting the differing characteristic helical geometry of dsDNA and dsRNA, whereas DRHs exhibit mixed behaviour depending on DNA-strand pyrimidine content. CTCF binding sites exhibit a distinct groove width signature. Persistence-length spectra from [~] 9 million sequences indicate that dsRNA is stiffer than dsDNA, whereas DRH exhibit intermediate stiffness modulated by DNA strand pyrimidine content. Persistence length increases upon CpG-modification, but decreases on hypermodification. Overall, the cgNA+ model enables a first, highly accurate, very large-scale, comparative study of sequence-dependent mechanics both within and across dsNA classes, demonstrating previously hidden regulatory layers. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=74 SRC="FIGDIR/small/740023v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@6de6acorg.highwire.dtl.DTLVardef@1435181org.highwire.dtl.DTLVardef@9c2b60org.highwire.dtl.DTLVardef@e3abf0_HPS_FORMAT_FIGEXP M_FIG C_FIG